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Abstract 

This is a continuation of two recent publications of the authors mi n] about 
reconstruction procedures for 3-d phaseless inverse scattering problems. The main 
novelty of this paper is that, unlike jT8], the Born approximation for the case of 
the wave-like equation is not considered. It is shown here that the phaseless inverse 
scattering problem for the 3-d wave-like equation in the frequency domain leads to 
the well known Inverse Kinematic Problem. Uniqueness theorem follows. Still, since 
the Inverse Kinematic Problem is very hard to solve, a linearization is applied. More 
precisely, geodesic lines are replaced with straight lines. As a result, an approximate 
explicit reconstruction formula is obtained via the inverse Radon transform. The 
second reconstruction method is via solving a problem of the integral geometry 
using integral equations of the Abel type. 

Keywords: phaseless inverse scattering, wave equation, reconstruction formula, Radon 
transform 

AMS classification codes: 35R30. 


1 Introduction 

The Phaseless Inverse Scattering Problems (PISPs) arise in applications to imaging of mi¬ 
crostructures of sizes of the micron range of less (1 micron=l/zm = 10 _6 m). In particular, 
this includes imaging of nano structures of sizes of hundreds of nanometers (~ 10 -7 m) 
and biological cells whose sizes are in the range of (5,100) /im [221 30] - To image these 
objects, one needs to use either optical radiation with the wavelength less than lfim or 
the X-ray radiation. However, it is well known that only the intensity (i.e. the square 
modulus) of the corresponding complex valued wave field can be measured for such small 
wavelengths. The phase cannot be measured, see, e.g. [sununiEsiES]. Therefore, we 
arrive at the problem of the reconstruction of the scatterer using only the intensity of the 
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scattered wave field. In this case the propagation of the wave held is governed by the 
wave-like PDE in the frequency domain. 

A similar problem, although for the Schrodinger equation in the frequency domain, 
arises in the quantum inverse scattering, where only the differential scattering cross-section 
can be measured, which is actually the square modulus of the solution of that equation, 
see, e.g. page 8 of [23] and Chapter 10 of [5j. Unlike the wave-like equation (12.6p . in the 
case of the Schrodinger equation the function f3 (x) (see (12.61) 1 is not multiplied by k 2 , 
which simplifies the problem. Note that, unlike PISPs, the conventional inverse scattering 
theory is based on the assumption that both the intensity and the phase of the complex 
valued wave field are measured, see, e.g. j3[ j8] [23, 24] . 

A reconstruction procedure for a 3-d PISP for the wave-like PDE in the frequency 
domain was proposed by the authors in [18]. In [T8] the linearization based on the Born 
approximation was used. However, the Born approximation assumption is inconvenient, 
since it is actually assumed there that k 2 \ j3 (a;)| << 1, where f3 (x) is the unknown coef¬ 
ficient and k is the frequency. Hence, the Born approximation breaks down for k » 1. 
On the other hand, large values of k are used in the reconstruction formula of [18]. Thus, 
the goal of this paper is to lift the assumption about the Born approximation. We still 
linearize the problem. However, the main novelty here is that our linearization does not 
break down when k —> oo. We achieve this via an extensive use of the structure of the 
fundamental solution of an associated hyperbolic equation with a variable coefficient in 
its principal part. The latter is the main new element here as compared with mm- 
This structure was derived in [33], also see [32]. 

The first uniqueness result for the PISP was proven in [13] for the 1-d case with some 
follow up publications [BCE]. Uniqueness for the 3-d case was proven in m m m- 
However, proofs in these references are not constructive. In fact, besides just uniqueness 
only it is desirable to develop rigorous numerical methods for the phase reconstruction. 
Many heuristic attempts were made by physicists to reconstruct the phase, see, e.g. [5J 

El EH m, 135]. 

However, rigorous numerical methods for PISPs were derived only very recently by 
the authors and Novikov [1?, [18| [26, 27] . In [17] a 3-d PISP for the Schrodinger equation 
was considered and an explicit reconstruction formula was derived, which is based on the 
inverse Radon transform, see, e.g. [21] for this transform. Thus, a long standing problem 
posed in 1977 in Chapter 10 of [3] was addressed in mi for the first time. In [T8J the 
result of ra was extended to the case of the wave-like equation (12.61) using the Born 
approximation assumption. 

While only unknown scatterers are involved in measurements in p~Tl ITS] , the recon¬ 
struction formula of [26] requires the involvement of two more known scatterers. A quite 
general reconstruction procedure of [27] is using measurements of the intensity of the full 
wave field on at least two spheres in the far field zone. The latter is unlike [T73> CEB) 26], 
where the intensity of only the scattered wave field is measured on just one surface and 
the far field approximation is not used. Note that usually the intensity of the scattered 
rather than the full wave held is measured. This is because measurements are conducted 
only outside of the brightening area since detectors are “burned” in the brightening area. 
For example, in images on page 22 of [6] the areas with brightening are depicted in the 
red color. On the other hand, outside of the brightening area the intensity of the full 
wave held is well approximated by the intensity of the scattered wave held. While works 
nn na Eg El EBl IZS EZ| are concerned with the reconstructions of coefficients of PDEs 
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from the phaseless scattered data, publications PHD] consider numerical methods of the 
reconstruction of shapes of obstacles from the phaseless scattered data. 

We present here two reconstruction methods for two PISPs. In both cases we use 
measurements of the intensity on a single sphere on many frequencies and the far field 
approximation is not used. First, we reduce our PISP to the Inverse Kinematic Problem 
IU1EDE2]. This leads to a new uniqueness result, which is based on two items: (1) a 
uniqueness theorem of [32] for the Inverse Kinematic Problem and (2) our reconstruction 
procedure. Next, we linearize the Inverse Kinematic Problem, as in [201 311 [32j . and 
reconstruct the unknown coefficient via the inverse Radon transform, as in mm - in 
the second approach we obtain after that linearization a problem of the integral geometry 
and solve it explicitly via solving certain integral equations of the Abel type. In our 
linearization we assume, similarly with [2Di i3T] [32], that certain integrals over geodesic 
lines are actually integrals over straight lines. 

As it is often happen to other reconstruction procedures (see, e.g. [26] EZ]) , our two 
reconstruction procedures cannot be considered as algorithms yet. In other words, they 
cannot be considered as sequences of steps leading to the numerical solutions. To turn 
them in algorithms, our reconstruction steps must be complemented by some regulariza¬ 
tion steps. Although the latter is possible, we leave this to future numerical publications 
for brevity. 

In section 2 statements of problems under consideration are presented. In section 3 
we consider an auxiliary Cauchy problem for a hyperbolic equation, which is important 
for our study. In section 4 we present a reconstruction method via the inversion of the 
Radon transform. Finally in section 5 we present the second reconstruction method via 
solution of a problem of the integral geometry. 

2 Problem statement 

Let B > 0 be a number. Let fl = {|x| < A} C M 3 be the ball of the radius R < B 

with the center at {0} Denote Y — {|x| < B} , S — {|x| = B} . Let n(x),x 6 I 3 be a real 

valued function satisfying the following conditions 

n G C 15 (M 3 ), 7i 2 (x) = 1 + P(x), (2.1) 

P(x) > 0, P(x) = 0 for x 6 l 3 \ fl. (2.2) 

The smoothness requirement imposed on the function n(x) is clarified in the proof of 
Theorem 1 in subsection 3.1. The conformal Riemannian metric generated by n(x) is 

dr = n(x) \dx \, \dx\ = \JJdx i) 2 + (dx 2 ) 2 + (dxs) 2 . (2.3) 

Below we impose the following Assumption: 

Assumption. We assume that geodesic lines of the metric (j2.3j) satisfy the regularity 
condition, i.e. for each two points x, y € M 3 there exists a single geodesic line F (x, y) 
connecting these points. 

It is well known from the Hadamard-Cartan theorem [2] that in any simply connected 
complete manifold with a non positive curvature each two points can be connected by a 
single geodesic line. The manifold (fl, n) is called the manifold of a non positive curvature, 
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if the section curvatures K(x,a) < 0 for all x G Q and for all two-dimensional planes a. 
A sufficient condition for K(x,a) < 0 was derived in 


<9 2 In nix) s 


i,j =1 


dxjdxj 


For x,y E R 3 , let r(x,y) be the solution to the problem 

|Vxr(a;,i /)| 2 = n 2 (x), t(x, y) = O (|x - y \), as y ->• x. (2.4) 

Let da be the euclidean arc length of the geodesic line V (x, y ). Then 

r{x,y)= I n{C)da. (2.5) 

r(x,2/) 

Hence, r(x,y) is the travel time between points x and y due to the Riemannian metric 
(12.31) . Due the Assumption, r(x,y) is a single-valued function of both points x and y in 
R 3 x R 3 . 

We consider the following equation 

Am + k 2 n 2 (x)u = — 5(x — y), x G R 3 , (2.6) 


where the Laplace operator is taken with respect to x, the frequency k > 0 is a positive 
real number and y G R 3 is the source position. Naturally, we assume that the function 
u(x, k,y) satisfies the radiation condition 

Q u 

——K iku = o(r _1 ) as r = lx — y\ —» oo. (2.7) 

or 

Denote u 0 (x, k,y ) the solution of the problem ( 12 . 6 p . (12.7[) for the case n(x) = 1. Then u 0 
is the incident spherical wave, 


uq(x, k, y) 


exp (—ik \x — y\) 
47t \x — y | 


Let u sc (x, k, y) be the scattered wave, which is due to the presence of scatterers, in which 
n(x) 7 ^ 1. Then 


exp (— ik \x — y\) 

u sc \x, k, y) = u(x, k, y) - u 0 {x, k, y) = u(x, k, y) - — -:-. ( 2 . 8 ) 

47t \x — y | 

Combining Theorem 8.7 of (4j with Theorem 6.17 of [7] and taking into account that n G 
C 15 (R 3 ), we obtain that the problem (12.61) . (12.71) has unique solution u G C l&+a (\x — ar 0 | > rj) , Vr/ > 
0 for any a G (0,1). Here C k+a are Holder spaces. 

We consider the following two Phaseless Inverse Scattering Problems (PISPs): 

PISP1. Suppose that the following function f\ (x, k, y) is known 

fi {x, k, y) = | u sc (x, k, y) | 2 , V (x, y) G S x S, Wk > k 0 , (2.9) 


where k 0 = const. > 0. Determine the function (3 (x) satisfying conditions 112. 1\) . H2.2 1) . 


4 












PISP2. Suppose that the following function f 2 (x, k, y ) is known 

f 2 {x, k, y ) = | u(x, k, y) | 2 , V (x, y) E S x S, Vk > k 0 , (2.10) 

where k 0 = const. > 0. Determine the function f3 (x) satisfying conditions h2.1\) . \2.21) . 

In both these problems we do not take into account the fact that the intensity of the 
wave field can be measured only outside of the brightening area (see Introduction). This 
case is a more difficult one than we now consider and we hope to study it in the future. 

Below we present a reconstruction method for the PISP1 in a linearized case. This 
method is based on the inverse Radon transform. In addition, we present a different 
reconstruction procedure both for the PISP1 and for the PISP2 (both are linearized) via 
solving a problem of the integral geometry. Note that while the inverse Radon transform 
is applicable to the PISP1, it is inapplicable to the PISP2. On the other hand, the 
reconstruction methods based on the integral geometry work for both these problems. 

Remark. In fact, it follows from our reconstruction procedures that, in the lin¬ 
earized case, it is sufficient to know functions fi (x, k, y ), f 2 (x, k, y) only for all (x, y) E 
(S fl |:r 3 = z}) x (S' 0 {x 3 = z}) for each 0 E (-R, R) and for all k > k 0 . This case of 
data collection is more economical than the one of (12.91) . (12.101) . In this case the function 
fd ( x ) is reconstructed separately in each 2-d cross-section hi fl (x 3 = z } of the domain 12. 


3 The auxiliary Cauchy problem for a hyperbolic equa¬ 
tion 

Consider the following Cauchy problem 

n 2 (x)v t t = An + 8{x — y, t), xEM. 3 ,t> 0, (3.1) 

v (x, 0) = v t (x, 0) = 0. (3.2) 


3.1 The form of the solution of the problem (13.11) . (13.211 


Let ( = (Ci 5 C 25 C 3 ) 5 C — C( x > y) be geodesic coordinates of a variable point x with respect 
to a fixed point y in the above Riemannian metric (12.31) . By the above Assumption, there 
exists a one-to-one correspondence x <=>• ( for any fixed y. Therefore, for any fixed point 
y the function ( = £( x,y) has the inverse function x = f (C,y) which determines the 
geodesic line T(x, y) = {£ : £ = / (s( 0 , y) ,s E [0, r(x, ?/)]}, where ( 0 is the vector which is 
tangent to T (x,y) at the point y, is directed towards the point x and also |£ 0 | = n~ 1 (y). 
Moreover, by the formula (2.2.28) in the book [32j the function ((x,y) can be expressed 
via the function r(x,y) as 


C(*. y ) = v v t2 ( x ’ y)- ( 3 - 3 ) 

Note that in our case one should take in the formula (2.2.28) of [23] A(y) = n' 2 (y)J, 
where / is the unit matrix. Also, since by (12.11) n(x) E C 15 (R 3 ) and the Assumption 
holds, then T 2 (x,y) E C 15 (M 3 x M 3 ), k > 2, and ((x,y) E C 14 (M 3 x M 3 ) [33]. Consider the 
Jacobian J(x, y), 

J(x,y) = det|^. (3.4) 
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By the formula (2.2.18) in [33J ( = x — y + 0(\x — y | 2 ) as x —> y. Hence, 


d(i ,. hS^ + Oih 2 ) 

7— = lim - - - 

UXj x= y x j Vj—h^O h 


^iji hO — 1)2, 3, 


where Sij is the Kronecker delta. Hence, J(y, y) = 1. Next, the smoothness of the function 
C(x,y) with respect to x implies that the function J(x,y ) is continuous with respect to 
x. Further, since the correspondence x ( is one-to-one and differentiable, then (13.411 
implies that J(x,y ) 7 ^ 0, V (x,y) 6 l 3 x M 3 . Hence, 


J(x, y) > 0, V ( x , y) e M 3 x M 3 . 


(3.5) 


For any y G M 3 and for any T > 0 consider two domains D(y, T) and D*(y , T), 


D (y, T ) = {(x,t) : 0 <t< T-\x-y\} and D*(y,T) = {(x,t) : r(x,y) < t < T-r{x,y)}. 

Theorem 1. Let the function n(x) satisfies conditions H 2. 1\) and \2.2 1) and let the 
Assumption holds. Then for any T > 0 and for any fixed y G M 3 there exists unique 
solution v(x,t,y ) to the problem \ 3. 1\) . \3.2\) . which can be represented in the domain 
D(y,T) in the form 

v(x, t, y) = A{x, y)8(t - r(x, y)) + v{x, t, y)H(t - r(x, y)), (3.6) 


where A(x,y) > 0 is given by the formula 

n 2 (y)V J ( x ’y) 


A{x,y) = 


47 m(x)r(x, y) 


(3.7) 


H(t) is the Heaviside function and the function v(x,t,y) G C 2 (. D*(y,T )). 

Proof. We note first that we can use a/ J(x, y) in (13.711 since by (I3.5|l J(x,y) > 0. 
Rewrite the problem (13.ip . (j3.2j) in the form 

v tt — div (n" 2 (x)Vr) + Vn^ 2 (x) • Vn = n~ 2 {y)5{x — y,t), n| t<0 = 0, (3.8) 

We now use results obtained in the book [33]. Represent the solution to the problem 
for t > 0 in the form 

v (x, t, y) = 

S 

a- i(x, y)8{t 2 - t 2 (x, y)) + a k (x, y)0 k (t 2 - t 2 (x , y)) + w 8 (x, t, y) 


n 2 (y) 


k =0 


(3.9) 


where the integer s > 1 is chosen below and functions d k (t) are defined as 

Qk(t) = —H(t), k = 0 , 1 , 2 ,.. .,s. 

Coefficients a k (x,y) in (13. 9 p are dehned by the following formulas 
^/J(x,y)n 4 (y) 


a-i(x,y) = 
a r (x, y) = 


27 xn(x) 
a~i(x,y) 


(3.10) 


4 T r ~ l {x,y) 


T r {f,y){a-i{f,y)) n (0 A €ffi-i(^, y)da, r > 0, (3.11) 


r (x,y) 
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where £ is a variable point on this geodesic line T(x, y). We have derived formulae (I3.10p . 
(13.lip from formulas (2.2.41)-(2.2.44) of the book [33], in which the notation a r was used 
for a r and dr' = n(£)dcr. We took m — 1 in the latter formulae. 

The residual function w s (x,t,y) in (13.91) is the solution of the following problem 

d 2 w s - n~ 2 (x)Aw s = F s (x, t, y ), w\ t<0 = 0, (3.12) 

F s (x,t, y) = n~ 2 (x)(Aa s (x,y))6 s (t 2 - r 2 (x,y)). (3.13) 

Let k > 1 be a sufficiently large integer which is chosen below. If the function n(x) G 
C k ( R 3 ), then r 2 (x,y) G C k (R 3 x R 3 ) (see [33], p. 27). Therefore it follows from (13.101) 
and (13.111) that 

a_i(x, y ) G C k ~ 2 { R 3 x R 3 ), a r (x, y ) G C fc - 4 " 2r (M 3 x R 3 ), r = 0,1,..., s. 


Hence, the function F s (x,t,y ) G C l (D(y,T)), where l = min(m — 6 — 2s, s) and this 
function vanishes for t < r(x,y). Hence, w s (x,t,y ) = 0 for t < r(x,y). 

Using theorems 3.2, 4.1, Corollary 4.2 and energy estimates of Chapter 4 of [T9], one 
can easily prove that w s G {D(y, T)) , see, e.g. Theorem 2.2 of [T2j for a similar result. 
Choosing m = 15 and s = 3, we obtain l = 3. Thus, w 3 G C 2 (D(y, T)). Since by (13.71) and 
(13.101) A(x,y) = a-i(x,y)/(27i 2 (y)r(x,y)), then A(x,y) > 0 and A(x,y) G C 13 (R 3 x R 3 ). 
Moreover, setting 


v{x,t,y) = 


n 2 (y) 


^a r (x,y)- - +w 3 (x,t,y ) 


r =0 


T: 


, (x,y) G D*(T,y), 


we see that formula (13.91) coincides with formula (13.6p and also that the function v(x, t, y) G 
C 2 (D*(T, y)). □ 


3.2 Connection with the problem (12.61) . ( 12.71) 


Let $ C R 3 be an arbitrary bounded domain. Lemma 6 of Chapter 10 of the book m as 
well as Remark 3 after that lemma guarantee that functions d k v(x, t, u), k = 0,1,2 and 
A X v(x,t,v) decay exponentially as t —> 00 , as longs as x reminds in the domain K. In 
other words, there exist constants M = M ($, j3) > 0, c = c (<E>, j3) > 0 such that 


\d k v (x,t,y )| , \Av(x,t,y)\ < Me ct for all t > t 0 and for all x G $, (3.14) 


where to = to (4*, f$) = const. > 0. Hence, one can consider the Fourier transform 
V (x, k, y) of the function v, 


OO 


V(x, k, y) 


J v (x, t, y) exp (— ikt) dt. 
0 


(3.15) 


Next, Theorem 3.3 of [[36] and Theorem 6 of Chapter 9 of E3 guarantee that V(x, k, y) = 
u(x, k,y), where the function u(x, k,y ) is the above solution of the problem (12.61) . (12.71) . 

Comparing (13.14p and (13.151) with (13.61) . integrating by parts in (13.151) and also using 
(12.81) . we obtain 


u (x, k, y) 


exp(— ikr(x, y)) 


A{x, y) + O 



, k —> 00 , 


(3.16) 
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,k—>oo. (3.17) 


/ , \ , / s , , .. exp (—ik \x — y\) /I 

u sc {x, k, y) = A{x , y) exp (-ikr(x, y)) - — ^ ^ -b O I - 


Hence, (12.9]) , (12.1 OH . (13. 161) and (13. 1 71) imply that for (x,y) e S x S 

fi (x,k,y) = 


A 2 (x,y) + 


1 _ A (x,y) 

167T 2 \x — y\ 2 27t \x — y\ 


cos [k (t(x, y) — \x 


+ O ( — ) , k —> oo, (3.18) 


h ip, k , v) = A 2 (x, y) + 0 ( - ) , k oo. 


(3.19) 


4 Approximate solution of the linearized PISP1 via 
the inverse Radon transform 

4.1 Reconstruction of the function r(x,y) for (x,y) E S x S 

Below in this subsection k > ki > /c 0 , where k\ » 1 is a sufficiently large number. We 
now fix the point (x,y) £ S x S and consider fi(x, k,y ) as the function of k for k > k\. 
ft is possible to hgure out whether or not r (x, y) = \x — y \. Indeed, it follows from (13.191) 
that r(x,y) = \x — y\ if and only if lim^oo fi(x,k,y) exists. Ignoring in (13.18)) the term 
O ( k~ l ), we obtain for k > k\ 


fi (x, k, y) = A 2 (x,y) + 


A(x, y) 


cos [k (t(x, y) — \x 


167T 2 \x - y\ 2 27t \x — y\ 

In (14.ip we use instead of . By (14.Tj) there exists a number k 2 > k { such that 

1 


(4.1) 


fi(x, y) = fi (x, k 2 , y) = max /, (x, k,y) = A(x, y) + . , , 

k>k i y 47r|rc — y\ 

Hence, we find the number A{x, y) as 

A(x,y) = \/f;(x,y) 


(4.2) 


47r|a; — y\ 


(4.3) 


Assume that r(x,y) ^ \x — y\. Then the positivity of the function f3 (x) and (12.5p 
imply that r(x,y) > \x — y\. Choose the number k 3 > k 2 such that 


k 3 = min {k : k > k 2 , f x (x, k, y ) = /*(x, y)} . 
Then (14.II) - (14.41) imply that 

h (t(x, y)-\x- y\) = k 2 (r(x, y) - \x - y\) + 2vr. 
Therefore, we reconstruct the number r(x, y) as 

2vr 


(4.4) 


x(x,y) = \x-y\ + 


k 3 -k 2 


(4.5) 


























Next, we should reconstruct the function (3 (x), which is done below. To do this, we 
consider a linearization. However, we now can formulate uniqueness theorem for PISP1 
without the linearization assumption. This is because the knowledge of the function 
r(x, y) for all (x, y) G S x S, which follows from (14.5)) . is equivalent to the so-called Inverse 
Kinematic Problem. This problem was studied in [20J [3TJ [32]. In particular, Theorem 3.4 
of Chapter 3 of [35] claims uniqueness of the reconstruction of the function n (x) from 
the knowledge of the function r(x, y) for all (x, y) G S x S. Thus, our Theorem 2 follows 
immediately from Theorem 3.4 of Chapter 3 of [32j as well as from (14.51) . 

Theorem 2 (uniqueness). Suppose that two functions n\ (x), ri 2 (x) satisfying condi¬ 
tions H2.1\ ) . \2.2 1) also satisfy Assumption of section 2. In addition, assume that these two 
functions generate the same function f\ (x, k , y) in H2.9\) and that the function fi (x, k , y) 
has the form i.e. the term O (1/k) in \3.18\) is dropped. Then ni (x) = n 2 (x). 

4.2 Reconstruction of the function (3 (x) 

Assume that 

ll^llc 2 (n) << 1- (4-6) 

Then the linearization of the function r (x, y) with respect to the function 3 leads to 

r (x, y) — \x — y \ + J 3 (f)da, (4.7) 

L(x,y) 

where L (x, y) is the segment of the straight line connecting points x and y and da is its 
arc length, see Theorem 11 in Chapter 3 of EDI as well as §5 in Chapter 2 of EH and §4 
in Chapter 3 of [32]. To be precise, in (14.71) one should have instead of “=". Since 
the function r (x, y) was approximately reconstructed via (j4.5[) . then we obtain from (14.71) 
that the following function h (x, y) — r (x, y) — \x — y\ is known 

h (x, y) — f P(Qds, V(x,j/) G S x S. (4.8) 

L{x,y) 

For any number z G R consider the plane P z = {X 3 = z }. Consider the disk Q z = 
Y A P z and let S z — S D P z be its boundary. Clearly Q z 7 ^ 0 for z G ( —B , B) and Q z = 0 
for \z\ > B. Denote 0 Z = (0,0, z) G Q z the orthogonal projection of the origin on the 
plane P a . For an arbitrary z G (—B,B) denote B z = V B 2 — z 2 the radius of the circle 
S z . We have 

B B 

Y= |J Q z , dY S — (J S z ,QcY. 

z=—B z=—B 

We now introduce some notations of the Radon transform, which we take from the 
book [5T]. Since our reconstruction formula is based on the inversion of the two-dimensional 
Radon transform, we now parametrize L (x, y) in the conventional way of the parametriza- 
tion of the Radon transform Ezn- Let v be the unit normal vector to L (x, y) lying in the 
plane P z and pointing outside of the point 0 2 . Let a G (0, 27t] be the angle between v and 
the Xi—axis. Then v = v (a) = (cos a, sinct) (it is convenient here to discount the third 
coordinate of v, which is zero). Let s be the signed distance between L(x,y) and the 
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point 0 2 (page 9 of [21]). It is clear that there exists a one-to-one correspondence between 
pairs (x, y ) and (v (a ), s), 

(x, y ) (v (a), s ), (x, y) G S z x S z ,a = a (x, y) G (0, 27t] , s = s (x, y) G (~B Z , B z ). 

(4.9) 

Hence, we can write 

L(x,y) = {r, a = (ri,r 2 ,a) : (r, z/(a)) = s}, (4.10) 

where r = (r'i, r 2 ) G M 2 , (,) is the scalar product in M 2 and parameters a = a (x,y) and 
s = s (x, y) are defined as in 04.91) . 

Consider an arbitrary function g = g (r) G C A ( P z ) such that g (r) = 0 for y G P z \Q z - 
Hence, 

J g(r) da = I g(r)da,V(x,y) G S z x S z , (4.11) 

L(x, y) (r,L'(a))=s 

where a = a (x, y ), s = s (x, y) are as in 04.9p . The parametrization of L ( x , y) in 04.lip 
is as in Q4.10p . Therefore, using 04.9p - 04.lip . we can define the 2-d Radon transform Rg 
of the function g as 

(Rg){x,y) = (Rg)(a,s) = J g(r)da,V(x,y) G S z x S z . 

(r,u(a))=s 

Hence, by 04. 8p and 04. 9 p 

h (x, y) = (Rfd) (a, s) = J (5 (r, a) da, V (x, y) E S z x S z ,Vz G (-B, B). (4.12) 

( r,v(a))=s 

Finally, since the function h ( x , y) is known, then 04.12p leads to the following reconstruc¬ 
tion formula 

P Or,r 2 , a) = i?” 1 (h(x,y)) (r 1 ,r 2 ,a ), (. x,y) G S z x G (~B,B). (4.13) 

Formula 04.13p is our final reconstruction result for PISP1 via the inverse Radon transform. 
Since the formula for R” 1 is well known [21], we do not specify it here for brevity. 

We now formulate uniqueness result for the linearized PISP1. This result follows 
immediately from 04.12p . 

Theorem 3 (uniqueness). Suppose that two functions (x ), ft 2 i x ) satisfying con¬ 
ditions \2.2 j) . Also, assume that corresponding functions n\ ( x ) ,n 2 ( x ) satisfy As¬ 

sumption of section 2. In addition, assume that these two functions generate the same 
function fi ( x , k , y) in (fiOjl and that the function fi ( x , k, y) has the form m i.e. the 
term 0(1 Ik) in is ianored. Finallu. assume that the linearization (O is valid. 

Then (x) = fi 2 ( x )• 

5 Reconstruction via the integral geometry 

In this section we present the reconstruction method for the function f3(x) from the 
function A(x,y ) given for all (x, y) G S x S. This method works for both PISP1 and 
PISP2. Indeed, consider two cases: 
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Case 1: PISP1. Given the function fi ( x , k, y) in and dropping the term O (1/k) 
in 03.181) . we obtain the function A (x,y) via 04.3p . 

Case 1: PISP1. Recall that by 02.11) and 02.21) n(x) — 1 for x £ S. Hence, in this 
case 02.101) . 03.71) and Q3.19|) imply that 

lim v 7 h Oh K V ) = A ( x , y) = Y J [ X ' V \ i V(x,y) e S x S, (5.1) 

k-* oo A7TT[X,y) 

where the determinant J(x,y ) was defined in 03.41) . 


5.1 Derivation of the problem of the integral geometry 

We again use 04.61) and the linearization 04.71) . Rewrite 04.71) as 


i 


r(x,y) = \x-y\ 


1 + J /3{y + s(x - y))ds 


(5.2) 


Because of 04.61) . we ignore in our approximate formulas below all terms containing the 
square of the integral in 05.21) as well as products of its derivatives. First, using 03.41) . we 
find an approximate formula for the determinant J(x,y) in 05.11) . Using 03.3p . we obtain 
the following approximate formula for ((x,y) 


C{x,y) 


{x - y) 
n 2 (y) 


1 + 2 / f3(y + s(x 



(5.3) 


x -^wS v ^ v+s(x 

o 


By (E3D we have for i,j = 1, 2, 3 


2 /))(l - s)ds. 


dCifay) 

dxj 


n 2 (y) 


1 + 2 J (3(y + s(x 

o 


y))ds 


i i 

(Xi-Vi) J P Xj {y + - y))sds - (Xj -Vj) J p x .(y + s(x - y))(l - s)ds 

0 0 
^ 1 

- ^ 2 J) / P XiXj (y + s(x~y))s(l-s)ds. 

0 

Hence, in the 3x3 determinant J(x,y ) = det (d(/dx) the product of diagonal terms 
dominates the rest of terms, which should be set to zero because of products of the above 
mentioned integrals with the function ft and its derivatives. Hence, an approximate 
formula for J(x, y) is 


2 

n 2 (y) 


J(x,y ) 


n 6 (y) 


1 + 6 / P(y + s(x 


y))ds + 2(x 


y) J WP(y + s(x-y))(2s-l)ds 
o 
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(5.4) 


\x — y | 2 / A /3(y + s(x — y))s( 1 — s)ds 


Note that 


(x-y)VP(y + s(x-y)) = 


d(3(y + s(x — y)) 
ds 


Hence, 

1 i 

(x-y) [ V/% + s(x - y))(2s - 1 )ds = j3(x) + jd(y) - 2 [ P{y + s{x - y))ds. 


Hence, we obtain 


J(x,y) = 1 + 2 / P(y + s(x- y))ds 


\x — y | 2 / A/3(y + s{x — ?/))s(l — s)ds , (x, y) E S x S. 


(5.5) 


Next, by fl5T 


r(x,?/) |x-y| 

Thus, fl5.ll) . (15.4p and (15.51) imply that 


1- / P{y + s{x - y))ds , 


(5.6) 


A(x,y) = 


47r|x — 7/1 


1 - 


\x - 7/p 


A f3(y + s(x — y))s( 1 — s)ds , (x,y) E S x S. (5.7) 


Both in (15.61) and (15.71) we again use ’ instead of Denote A /3(x) = q{x). Then 
we obtain the following problem: 

The Problem of the Integral Geometry. Find the function q E C (Y) , q (x) = 0 

for x E T\0, assuming that the following integrals g(x,y) are given over segments 
L (x, y ) of straight lines connecting any two points (x,y) E S x S,x ^ y 


\x-y\ / q(y + s(x-y))s(l-s)ds = g(x,y), 


(5.8) 


g(x,y) =-8nA(x,y)+ 2\x-y\ \ 

Note that g(x, y) = g(y, x). Since the weight function s(l — s) is present in the integral 
(15.8p . the problem (15.81) cannot be solved via the inversion of the Radon transform as in 
section 4.2. Therefore, we propose a different method in subsection 5.2. If the problem 
(15.8p is solved, then the function (3(x) can be found via the solution of the Dirichlet 
boundary value problem for the Laplace equation, 

A /3{x) = q(x), x EY; /3(x)\ s = 0. 
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5.2 Solution of the problem of the integral geometry 


We now show that the function q(x) can be reconstructed from (15.8|) separately in each 
2-d cross-section Q z = Y Pi P Zl z G (-R, R) of the ball Y. First, we rewrite equation (15.8)1 
as 

\x-y\ 

/ ( (x — y \ \ 

q \ y + Sl \x - y\ ) Sl(y \x ~ y \ ~ s ^ dsi = 9(x,y), ( 5 . 9 ) 

o 

V(x,p) G S z x S z ,x ± y, 

where Si is the arc length of L (x, y). Recall that the plane P z = {x 3 = z}. We consider 
z € (—R,R), since q (x) — 0 for |x| > R (see (12.211 ). In the plane P z we introduce 
polar coordinates r, <p of the variable point £ = (£i,£ 2 ) as £ x = r cos<p,£ 2 = r sin <p. We 
characterize the segment of the straight line L(x, y) passing through points x,y G S z by the 
polar coordinates (p, a) of its middle point (x + y)/ 2. Hence, \x — y\ — 2 \JB 2 — z 2 — p 2 . 
Change variables in the integral (15.9)1 as 

Si r, Si = 1 / B 2 — z 2 — p 2 — \/r 2 — p 2 , 

Then si(|x — y\ — s 1 ) = B 2 — z 2 — r 2 . The equation of L(x, y) can be rewritten as 

(p = a + (—l)- 7 arccos —, j = 1,2, r > p, (5.10) 

where j = 1 corresponds to the part of the straight line L(x, y) between points (x + y)/2 
and x, and j = 2 to the rest of L(x, y). 

Note that ds\ = —rdr/\Jr 2 — p 2 . Obviously there exists a one-to-one correspondence, 
up to the symmetry mapping (x, y) ^ (y, x) between pairs (x, y) G S z x S z and (p, a) G 
(0, R) x (0,27 t). Denote q(r cos <p, r sin<p, z) = q(r,<p,z ) and g(x,y ) = g(p,a,z). Using 
(15.1011 . we rewrite equation (15. 9 p as 


2 

\ f /( r - a - H) 
j=i i 


s p \ ( B 

J arccos -, z) — 
r 


2 -z 2 


r 2 )rdr 


\Jr 2 — p 1 


= g(p,a,z), 


(5.11) 


where p 0 = y/R? — z 2 . Represent functions q(r, p, z) and g(p, a, z) via Fourier series, 


OO 

q(r,p,z)= q n (r, z) exp(in<p ), 

n =—00 


(5.12) 


OO 

g(p,a,z)= Y 9 n{p,z)exp(ina). (5.13) 

n =— 00 

Multiplying both sides of (15. lip by exp(— ina)/(2n) and integrating with respect to ct, we 
obtain for all n = 0, ±1, ±2,... 


Po 

J q n (r,z ) cos 

p 


n arccos : 


P 


2 (B 2 — z 2 — r 2 )rdr 

\Jr 2 — p 2 


9n(p,z). 


(5.14) 
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Denote 


Pn(r, z ) = q n (r, z){B 2 - z 2 - r 2 )r. 
Then equation (j5. 14j) becomes 
Po 


p n {r , zj cos l n arccos - 


p\ 2 dr 


rJ yV 2 - p 2 


= 9n(p,z), p G (0, p 0 \. (5.15) 


This is the integral equation of the Abel type. To solve equation (j5.15j) . we apply first 
the operator L to both sides of (15.151) . where 


Po 


L {h{p)) ( s ) = - 

7r 


1 f h{p)pdp 


\Tp 


t 2_ s 2 


,s G (0,p 0 ). 


Then changing the limits of the integration, we obtain 
Po 


~ / Pn(r,z) 


7T 


2 P ( P\j 

— . - . cos n arccos - dp 

V rJ 


dr — L (g n (p, z)) (s ). (5.16) 


Change variables in the inner integral (I5.16p as 

p 6, p 2 = s 2 cos 2 {0/2) + r 2 sin 2 {0/2). 


Then 

2 pdp = (r 2 — s 2 ) sin 0 cos OdO , 

\Jp 2 — s 2 • y/ r 2 — p 2 = {r 2 — s 2 ) sin 0 cos 0. 
Hence, equation (15.161) can be rewritten as 

Po 

/ p n {r , z)Q„(r, s)dr = L ( g n {p , z)) (s), 


(5.17) 


^ s 1 } ( V 7 " 2 cos 2 (0/2) + s 2 sin 2 (0/2) ^ 

*1 V r ) 


We have Q n {s,s) = 1. Hence, differentiating (15.17[) with respect to s, we obtain Volterra 
integral equation of the second kind 


Po 


d 


Pn{s 1 z )~ Pn{r , z)T n {r , s)dr = — — [L (^ n (p, 2)) (s)] , s G (0, p 0 ) , 


(5.18) 


T n (r,s) = 


ns 


Ti^/r 2 — 


sm n arccos 


y/r 2 cos 2 {0/2) + s 2 sin 2 (0/2) 


(5.19) 


x 


sin(6 , /2)c76 l 


y/r 2 cos 2 (0/2) + s 2 sin 2 (0/2) 
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It follows from (15.191) that the kernel of the Volterra integral equation (15. 18j) has the form 


r, s = 


T n (r,s ) 
\/r 2 - s 2 


where the function T n (r, s) is continuous for 0 < s < r < p 0 . Therefore, it follows 
from the theory of Volterra integral equations of the second kind that for each z £ 
(—R,R) there exists a solution p n (s,z ) £ C[0,p 0 ] of equation (j5.18j) and this solution 
is unique. Furthermore, it is well known from that theory that equation (15.18j) can be 
solved iteratively as 


Pn{ S , Z ) = 


d_ 

ds 


Po 


Pn(s,z) = 


Pt 1 


[L (. g n (p,z )) (s)], 


d 


(r, z)T n (r , s)dr - — [L (g n (p, z)) (s)] , k = 1, 2,... 


(5.20) 


(5.21) 


and this process converges in the space C [0, p 0 ] to the solution p n (s, z) of equation (j5.18j) . 
Formulae (I5.20|) and (15.211) finish our second reconstruction procedure. 

As a corollary, we formulate the following uniqueness theorem, which follows immedi¬ 
ately from the above reconstruction process. 

Theorem 4 (uniqueness). Assume that the function A (x , y) has its approximate form 
& Suppose that this function is given for all (x , y) £ S x S. Then there exists at most 
one function f3 £ C 2 (V) — 0 in Y\Q which is involved in fl 5. r /\ ). In particular, 

PISP2 has at most one solution as long as the function A (x , y) is as in The same 

is true for PISP1, if, in addition, the term O (l/k) is dropped in KS.lh 1) . □ 
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